Two- Site Kondo Effect in Atomic Chains 
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Linear CoCunCo clusters on Cu(lll) are fabricated by means of atomic manipulation. They 
represent a two-site Kondo system with tunable interaction. Scanning tunneling spectroscopy reveals 
oscillations of the Kondo temperature Tk with the number n of Cu atoms for n > 3. Density 
functional calculations show that the Ruderman-Kittel-Kasuya-Yosida interaction mediated by the 
Cu chains causes the oscillations. Calculations find ferromagnetic and antiferromagnetic interaction 
for n = 1 and 2, respectively. Both interactions lead to a decrease of Tk as experimentally observed. 
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Magnetic atoms with partially filled d or f shells in- 
duce strong electron correlations, which cause spectac- 
ular effects such as Mott metal-insulator phase transi- 
tions, heavy- fermion behavior or the occurrence of high 
temperature superconductivity. The rich physics of these 
systems is due to the interplay of local and nonlocal cor- 
relation effects. Local correlations are due to Coulomb 
interaction which makes the probability of an electron to 
hop into an unoccupied dor f orbital at a site i depend on 
the number of electrons already present at that site. Non- 
local correlations arise from electrons propagating from 
one site i to another site j with Coulomb scattering at 
both sites. In this way, the Ruderman-Kittel-Kasuya- 
Yosida (RKKY) interaction [1] can cause magnetic cor- 
relations between distant atoms and its interplay with 
(local) Kondo physics is believed to be responsible for 
magnetic instabilities in heavy- fermion systems. 

Two interacting magnetic atoms are the smallest pos- 
sible solid state structure to display nonlocal magnetic 
phenomena due to these correlations. They are conse- 
quently an important model system and may also be 
viewed as a keystone in a bottom-up approach to mag- 
netic nanotechnology. In single-impurity systems, cor- 
relations between delocalized band electrons and more 
localized electrons at the impurity define a characteristic 
low-energy scale referred to as Kondo temperature Tk, 
and manifest themselves as an Abrikosov-Suhl resonance 
of width keTx near the Fermi level. Interactions between 
two impurities may alter this energy scale and in the limit 
of strong interactions lead to new regimes including fer- 
romagnetically locked impurity spins or an interimpurity 
singlet [2l,|3|. 

In this work the characteristic low-energy scale is ex- 



plored as a function of the interaction between two impu- 
rities. A Co atom on a Cu(lll) surface is connected to a 
second Co atom by a chain of Cu atoms. By evaluating 
the widths of the Abrikosov-Suhl resonances from scan- 
ning tunneling spectra we find that Tk is significantly 
lower for short CoCu^Co (n = 1, 2) than for longer chains 
(3 < n < 6). Starting from n = 3, Tk exhibits even-odd 
oscillations with the number of Cu atoms. Combining 
density functional theory (DFT) and perturbative renor- 
malization group (RG) calculations we show that the ob- 
served behavior is due to magnetic interactions between 
the Co adatoms. The observed oscillations of Tk are 
shown to result from the RKKY interaction mediated by 
the Cu chains. Overall, both ferromagnetic and antifer- 
romagnetic exchange interactions between the impurities 
are shown to reduce Tk compared to the single-impurity 
Kondo scale Tk- 

A few scanning tunneling microscope (STM) exper- 
iments previously addressed magnetic interactions be- 
tween individual adsorbed atoms (adatoms). The mag- 
netic hysteresis of two interacting Co atoms on a Pt(lll) 
surface was investigated as a function of the interatomic 
distance [4]. Kondo physics was not addressed in that 
work. The antiferromagnetic coupling between a Kondo 
Co atom and an unscreened Fe atom on CU2N surface 
induced a splitting of the Kondo resonance similar to 
Kondo splitting of single atom in external magnetic field 
[5] . One Co adatom separated by a few Cu layers from an 
entire magnetic Co film can be similarly understood in 
terms of single-impurity Kondo physics [6]. In both cases, 
none of the quantum phase transitions associated with 
the two-impurity Kondo problem [7] can arise. A disap- 
pearance of the Kondo effect or a reduction of Tk in small 
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FIG. 2. (Color online) Tk of Co atoms in CoCunCo chains as 
a function of n (squares). A dashed line indicates Tk = 110 K 
of a single Co adatom attached to the end of a Cu4 chain. 
The circles (red) indicate calculated sm^/e-impurity Kondo 
temperatures, which would be expected for CoCunCo chains 
in the absence of magnetic Co- Co interactions. The sinusoidal 
line is a fit to data in the RKKY interaction regime with 
calculated periodicity. 



FIG. 1. (Color online) Left: Pseudo-three-dimensional repre- 
sentations of constant-current STM images (38 A x 21 A) of 
linear CoCunCo clusters (n — 1 ... 6). Right: dJ/dl/ spectra 
recorded above the Co atoms of the chains. Green dots are 
experimental data, solid lines indicate Fano line shapes fit to 
the data. The dashed lines are Fano lines for Tk = 110 K, 
which corresponds to a Co adatom attached to a Cu4 chain. 

C02 and Ni2 clusters have been reported [8|,[9|- Moreover, 
the Kondo effect of two Ni atoms on Au(lll) was found 
to be fully developed as soon as the Ni atoms were not 
in nearest-neighbor positions [9]. Finally, a broadening 
of the Abrikosov-Suhl resonance has been extracted from 
tunneling spectra of Co atoms in next-nearest neighbor 
positions and used to determine their exchange interac- 
tion E'ex However, the method used in Ref. [10] to 
relate Tk and T^ex has been questioned [H]. While in 
these cases the interaction between magnetic impurities 
was short-ranged and mediated by the substrate or by 
direct contact of the magnetic atoms, an intermediate 
chain of nonmagnetic atoms may affect the interaction 
[l2| . The present work focuses on the impact of such a 
chain on the Kondo effect. For the first time the long- 
range interaction of two impurities is explored over the 
entire range from weakly coupled Kondo impurities to 
the regime where the interimpurity exchange coupling 
and Tk are of the same order of magnitude. 

The experiments were performed with a STM operated 
at 7K and 10~^ Pa. Cu(lll) surfaces and chemically 
etched W tips were prepared by Ar+ bombardment and 
annealing. Single Co atoms were deposited at ~ 10 K us- 
ing an electron beam evaporator while single Cu atoms 
were transferred from the tip to the sample [13]. The 
adatoms were chemically identified by the presence (Co) 
or absence (Cu) of the Abrikosov-Suhl resonance. Spec- 
troscopy of the differential conductance (di/dF) was per- 
formed by a lock-in technique using a modulation ampli- 
tude of 1 mVrms- Only tips which reproduced the known 



spectrum of a single Co adatom were used. 

Linear clusters of Co and Cu atoms were fabri- 

cated by manipulation with the STM tip by first assem- 
bling Cun chains along the close-packed [110] direction 
and then attaching Co atoms (Fig.[TJ left column). For 
n > 3, Co atoms are readily discriminated from Cu atoms 
by their larger apparent height. The right column of 
Fig. [1] shows d//dy spectra (dots) acquired above the Co 
atoms of each CoCunCo chain. Spectra acquired at the 
two Co atoms were virtually identical. The Abrikosov- 
Suhl resonance appears as an indentation of the d//dy 
signal around zero voltage and appreciably broadens from 
CoCuCo to C0CU3C0. To quantify its width, the spectra 
were fit by a Fano line [ITI, dJ/dF oc ((7 + e)V(l -|- e^) 
(black line). Fit parameters were the asymmetry param- 
eter q and e = (eV — eK)/(kBTK) {V'. sample voltage, 
eK^ resonance energy). While ~ 0.1 and eK ^ 2meV 
are rather independent of the number of Cu atoms, Tk 
exhibits a pronounced variation (Fig. [2]). Tk nearly dou- 
bles from n = 1 (Tk ^ 46K) to n = 2 (^ 79 K), in- 
creases further to ^ 108 K for C0CU3C0, and then oscil- 
lates (n = 4: ^ 93K, n = 5: ^ 110 K, n = 6: ^ 91 K). 
The maxima of the oscillation match the Kondo tem- 
perature of a Co atom at the end of CoCus and C0CU4 
chains (Tk = 110 K), which approximate a CoCuoo chain 
(dashed line in Fig.[2j). 

To interpret the evolution of Tk a first-principles de- 
scription of the CoCu^Co chains on Cu(lll) has been 
developed. DFT calculations were performed using the 
Vienna Ab Initio Simulation Package [18] with the imple- 
mented Projector Augmented Waves basis sets [l9, l20| 
and employing a generalized gradient approximation 
(GGA) [2l| to the exchange correlation potential. The 
surface was modeled by (3 x 9) Cu(lll) supercells with 
slab thicknesses of three to six layers and the CoCunCo 
chains oriented along the [110] direction. The Kondo 
physics observed in the experiments are governed by Co 
3(i electrons which are highly sensitive to their atomic 
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FIG. 3. (Color online) (a) Local den- 
sity of states (LDOS) of Cu 4s and 
4p orbitals in a CoCusCo chain and 
in an infinitely long Cu chain (Cuoo) 
on Cu(lll). Cu sites at which the 
LDOS was evaluated are indicated by 
dots, (b) LDOS of Cuoo on Cu(lll) 
versus energy and wave vector along 
the chain direction, (c) Cross sec- 
tion of (b) at ^ = OeV. (d) Relaxed 
CoCunCo (3 < n < 5) chains with 
calculated magnetization density for 
ferromagnetically aligned Co spins. 
Light (red) and dark (blue) colors in- 
dicate up and down orientations, re- 
spectively. 



environment [iGj. Therefore, both adatoms and first sur- 
face layer atoms were relaxed until the forces fell below 
0.02 eVA~^ for each atom. 

First, we show that 5m^/e-impurity physics cannot ex- 
plain the experimental variation of Tk. Using the Co 
hybridization functions obtained from DFT in a RG ap- 
proach we estimated the variation of Tk with cluster 
size due to single-impurity effects like variations in the 
local density of states at the Co site [22|. Extrapo- 
lating from Tk = 110 K for CoCua and C0CU4 to the 
CoGunCo chains yields variations of Tk (Fig. [21 circles) 
which are much smaller than experimentally observed 
(Fig. [21 squares) and do not even follow the trend of the 
experimental data. 

Due to the reduced coordination Cu atoms within the 
chain move by 0.16 to 0.23 A towards the Gu(lll) sur- 
face, which corresponds to 8 to 11% of the (111) in- 
terlayer spacing. Go atoms at the ends of the chain 
move by 0.31 to 0.33 A toward Gu(lll). These relax- 
ations have a significant influence on the magnetic cou- 
pling between the Go atoms. This becomes clear from 
calculations of total energies of a GoGuGo chain for fer- 
romagnetic (tt) and antiferromagnetic (t|) spin orien- 
tations with and without geometric relaxation. The re- 
sulting exchange interaction energy, T^ex = E^i — T^tt' 
is 14 ± 6meV for relaxed GoGuGo chains, while it is 
significantly smaller |T^ex| ^ 1 meV without relaxations. 
From a calculation within the local density approxima- 
tion (LDA), Ee^ ^ 40meV was reported for the unre- 
laxed chain [12]. For relaxed G0GU2C0 chains we find 
^ex = — 17 ± 3 meV, which is close to an earlier LDA re- 
sult [12] and our GGA result T^ex ^ — 20meV for the un- 
relaxed chain. The signs of T^ex indicate ferro- and anti- 
ferromagnetic coupling between the Go atoms in GoGuGo 
and G0GU2C0 chains, respectively. 

Despite the ferromagnetic Go-Go coupling in the 
GoGuGo chain an Abrikosov-Suhl resonance with Tk = 
46 K is observed in the experiments, which may ap- 
pear surprising. The interplay of interimpurity coupling. 



H = —J Si • S2, and Kondo screening has been theo- 
retically addressed in terms of two-site spin- 1/2 Kondo 
models [H, For ferromagnetic interimpurity exchange 
interaction, J > 0, dominating over the single-impurity 
Kondo energy scale, keTk ^ J, locking of the two spins 
to a total spin. Si -\- S2 = followed by a spin-1 Kondo 
effect has been predicted [2]. In this case, the Kondo tem- 
perature is reduced, Tk ^ kufl/J Q. With Tk = 46K 
and Tk = 110 K, J is estimated as 23meV. As J is re- 
lated to ^ex = 14 ± 6 meV by a factor of the order of 
one, the experimentally observed reduction of Tk is well 
in line with T^ex as obtained from our calculations (2^ . 
Consequently, theory and experiment consistently sug- 
gest that GoGuGo is in the crossover region between two 
independent and two ferromagnetically locked Kondo im- 
purities [2, 3], where a narrowed rather than a completely 
suppressed Abrikosov-Suhl resonance is found. 

In the case of G0GU2C0, where the Go atoms couple 
antiferromagnetically, our GGA calculations show that 
|T^ex| ^ 2kBTK. In a particle-hole-symmetric case, the 
two-site spin- 1/2 Kondo model exhibits a quantum crit- 
ical point at J ~ 2.2kBTK separating ground states 
with an antiferromagnetically locked interimpurity sin- 
glet from two Kondo-screened impurities [3|. In sys- 
tems without particle-hole symmetry, the quantum crit- 
ical point is replaced by a crossover region, where the 
spectral weight of the Abrikosov-Suhl resonance is contin- 



uously reduced and evolves into a pseudogap feature [25- 
^28^. In this crossover region, there are two energy scales, 
Tl < Th^ characterizing the spin- fluctuations and the 
quasiparticle excitation spectra of the system 25, 26, 28| . 
The lower scale, T^, gives rise to the sharpest and most 
pronounced feature in the spectral function at the Fermi 
level [isl and characterizes the onset of local Fermi liquid 
behavior [28|. Hence, ^bTl should appear as the width 
of the experimentally observed Abrikosov-Suhl resonance 
and we refer to Tl = Tk as the Kondo temperature of 
the two-impurity system. As Tl < Tk [^2^, H, HH, a 
narrowed Abrikosov-Suhl resonance as observed experi- 
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mentally (Figs.[Tl [2j) is well in line with this crossover 
regime. The results for n = 1, 2 prove that positive and 



m2i 



negative exchange interactions lead to Tk < Tk 

To understand the oscillations of Tk for n > 3 the 
electronic structure of an infinite Cu chain (Cuoo) on 
Cu(lll) is considered. Its structural relaxations (i.e. 
0.2 A downward relaxation of the chain atoms) and its 
local density of states (LDOS) are similar to those ob- 
served from CoCunCo for n > 3. Cu atoms that are 
no direct neighbors of Co atoms exhibit a similar LDOS 
as Cu atoms in Cuoo [Fig.[3fa)]. Therefore, the Fermi 
wave vector /cp, which determines LDOS oscillations of 
the chains and the RKKY interaction, has been calcu- 
lated from the band structure of Cuoo- Figure [H^b) shows 
the energy- and momentum-resolved LDOS of conduc- 
tion electrons of Cuoo- The band of unoccupied states 
starting at ^ = 1.2eVat/c = has been measured in 
Refs. l3, [m and does not contribute spectral weight at 
the Fermi energy {E^ = OeV). At [Fig.[3l^c)] spec- 
tral weight is suppressed for kp < 0.1 (27r/a) [a: Cu(lll) 
lattice constant, a = 2.57 A] and exhibits a maximum at 
~ 0.37(27r/a). Hence, Co-Co RKKY interactions and 
LDOS resonances at Ep are expected to oscillate with 
a wave vector 2 /cp ^ 0.74(27r/a) which by subtracting 
a reciprocal lattice vector is identical with —0.26 (27r/a) 
and corresponds to a direct space period of ~ 3.8 a. This 
period can be clearly seen in Fig.[3jd) as an oscillatory 
magnetization density along the Cu chains. The period 
expected for Tk, however, is different. In the limit of 
weak RKKY interaction, the correction to the Kondo 
temperature reads f^-T^^ ^eJ^B [22", [H, |30| . Given 
that E^^ oc sin^(2 kpna) the spatial periodicity is reduced 
to ^ 1.9 a, which corresponds well to the even-odd oscil- 
lations of Tk observed in the experiments. 

In summary, linking two Co atoms by a chain of Cu 
atoms nonlocal correlations between two Kondo impuri- 
ties have been probed. The interimpurity interaction is 
proven to quench Kondo temperatures in short clusters 
and leads to RKKY induced oscillations at larger chain 
lengths. A reduction of the Kondo temperature indepen- 
dent of the sign of the interimpurity exchange interaction 
is found. These effects observed from a two impurity sys- 
tem may find a counterpart in crystalline solids as indi- 
cated by model studies of the double Bethe lattice [32]. 
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